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Abstract 

In this paper, we shall establish the spreading speed and existence of traveling waves 
for a non-cooperative system arising from epidermal wound healing and characterize 
the spreading speed as the slowest speed of a family of non-constant traveling wave 
solutions. Our results on the spreading speed and traveling waves can also be applied 
JLh 1 to a large class of non-cooperative reaction-diffusion systems. 
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1 Introduction 



In this paper, we study the spreading speeds and traveling wave solutions of a 
non-cooperative reaction-diffusion systems arising from wound healing. Wound 
healing is complex and remains only partially understood, despite extensive 
research. Several reaction-diffusion models have been developed in Sherratt 
and Murray [23f24] . Dale, Maini, Sherratt [6] and others to understand the 
biological process of epidermal wound healing through mathematical analysis 
and numerical simulations. We refer to Murray [19] for more detailed discus- 
sions and further references. The models consist of two conservation equations, 
one for the epithelial cell density per unit area {u\{x, t)) and one for the con- 
centration of the mitosis-regulating chemical (u2(x,t)). There are two types 
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of the chemicals, one in which the chemical activates mitosis and the other in 
which it inhibits it. The following simplified model was proposed in [23fT9] for 
the activator 

-— = diAui + s(u 2 )ui(2 - Mi) - ui 

d % (i-i) 

— - = d 2 Au 2 + b{h{u x ) - u 2 ) 
at 

where b > 0, k G (0, ~), s(u 2 ) = ku 2 + 1 — k is the linearized function which 
reflects the chemical control of motosis. The chemical production by the func- 
tion h{u\) = "H, 1V2 ' , C G (0,1) reflects an appropriate cellular response to 
injury. The qualitative form of the solution of (II. ip in the linear phase is of 
a wave moving with constant shape and speed. Such a solution is amenable 
to analysis if we consider a cone dimensional geometry rather than the two 
dimensional radially symmetric geometry. Mathematically, we look for a trav- 
eling wave solution of the form ui(x,t) = u\(^) , u 2 (x , t) = u 2 (£,),£, = x + ct 
where c is the wave speed, positive since here we consider waves moving to the 
left. In Section HI we shall establish the existence of traveling waves as well as 
the results on the speed of propagation to ( 11. II) . In addition, we characterize 
the minimum speed as the slowest speed of a family of non-constant traveling 
wave solutions of (II. ip . 

Traveling wave solutions and spreading speeds for reaction-diffusion equations 
have been studied by a number of researchers. Fisher [10] studied the nonlinear 
parabolic equation 

w t = w xx +w(l - w). (1.2) 

for the spatial spread of an advantageous gene in a population and conjec- 
tured c* is the asymptotic speed of propagation of the advantageous gene. His 
results show that (ll.2p has a traveling wave solution of the form w(x + ct) 
if only if |c| > c* = 2. Kolmogorov, Petrowski, and Piscounov [H] proved 
the similar results with more general model. Those pioneering work along 
with the paper by Aronson and Weinberger [Tp] confirmed the conjecture 
of Fisher and established the speeding spreads for nonlinear parabolic equa- 
tions. Lui [18] established the theory of spreading speeds for cooperative re- 
cursion systems. In a series of papers, Weinberger, Lewis and Li [TB]ll6f31|l52^ 
studied spreading speeds and traveling waves for more general cooperative 
recursion systems, and in particular, for quite general cooperative reaction- 
diffusion systems by analyzing of traveling waves and the convergence of ini- 
tial data to wave solutions. However, mathematical challenges remain because 
many react ion- diffusion systems are not necessarily cooperative due to var- 
ious biological or physical constrains. Thieme [2B] showed that asymptotic 
spreading speed of integral equations with nonmonotone growth functions can 
still be obtained by constructing monotone functions. For a related nonmono- 
tone integro-difference equation, Hsu and Zhao [12] , Li, Lewis and Weinberger 
[T7] extended the theory of spreading speed and established the existence of 



travel wave solutions. The author and Castillo-Chavez [29] prove that a class 
of nonmonotone integro-difference systems have spreading speeds and trav- 
eling wave solutions. Such an extension is largely based on the construction 
of two monotone operators with appropriate properties and fixed point the- 
orems in Banach spaces. A similar method was also used in Ma [22] and the 
author [28] to prove the existence of traveling wave solutions of nonmonotone 
reaction-diffusion equations. Weinberger, Kawasaki and Shigesada [35] discuss 
the minimum spreading speeds for a partially-cooperative system describing 
the interaction between ungulates and grass. It is cooperative for small pop- 
ulation densities but not for large ones. By employing comparison methods 
[35] established the spreading speeds of propagation. In a recent paper [30] . 
we study traveling waves and spreading speeds of propagation for a class of 
non-cooperative react ion- diffusion systems and a slightly different model de- 
scribing the interaction between ungulates and grass. 



In this paper, we shall study the spreading speeds and existence of traveling 
waves for the non-cooperative system arising from epidermal wound healing 
( II. ip . The minimum speed c* can be characterized as the slowest speed of a 
family of non-constant traveling wave solutions. In other words, we shall show 
that for c > c* (II. ip always has a nonconstant traveling solutions of the form 
(u\(x + ct),U2(x + ct)) with (mi(— oo),M2(— oo)) = but bounded away from 
zero at +oo, and there is no such traveling solution when < c < c*. Our main 
results for the epidermal wound healing model are summarized in Theorem 14. II 
The results for general non-cooperative systems (12. 3p are included in Theorem 
12.21 In Section H] we shall verify the assumptions of Theorem 12.21 for (11.11) and 
apply the general results to (11.11) . 

In order to better understanding of the spreading speeds and traveling solu- 
tions for the epidermal wound healing model. The general results in this paper 
have some significant improvements of those results over [30]. For example, in 
this paper we make use of the comparison principle from Fife [9]. Another form 
of the comparison principle from [35] was used in [30]. As a result, the assump- 
tions (H1-H2) and the proofs in Section [5] are somewhat different from [30j . 
By a suitable modification of the functions, the conditions in this paper seem 



easier to verify. In addition, Theorem 13. II is more general than that in [30], for 
example, d\ > d 2 is imposed in [30]. To take d\ < d 2 into consideration, some 
assumptions and proofs are substantially modified. In particular, for this epi- 
dermal wound healing model, we show that the condition for Theorem I3.1( ii) 
can be satisfied if ui(x) ^ 0. Finally, verifications of lower and upper solutions 
for the equivalent integral equations are significantly simplified via a result in 
Ma [21]. In [28] and [30], a more direct but lengthy verification of the lower and 
upper solutions are given for scalar and n-dimensional systems respectively. 
We also omit some standard proofs such as continuity and compactness for 
the operator which can be found in previous papers. 



2 Preliminaries 

We begin with some notation. We shall use Rjk,^, f, f^,r,u,v to denote 
vectors in M> N or iV-vector valued functions , and x, y, £ the single variable 
in R. Let u = (ui),v = (i>j) G R , we write u > v if Ui > Vi for all i; and 
u ^> v if Mj > f i for all i We further define for any r = (r J ) >> 0, r G R w the 
.R^-interval 

[0, r] = {m : < u < r, u G M^} C R N 

and 

C r = {u= («i) : Ui G C(R,R),0 < «i(ar) <r { ,xeR, i = l,...,N} 7 

where C(M., M) is the set of all continuous functions from R to R. 

Consider the system of react ion- diffusion equations 

ut = Du xx + f(u),x G R, t>0. (2.3) 

with 

u(x,0) = u (x), xGR, (2.4) 

where u = (w»), -D = diag(<ii, <i 2 , ■••, rfj\r)> rfj > 0, i = 1, ..., iV 

/(«) = (/l(«), f2(u), ..., /jv(w)), 

Mo( 2; ) is a bounded uniformly continuous function on R. In this paper, by a 
solution we mean a twice continuously different iable function u(x,t) in R x 
(0, oo) and continuous in R x [0,oo), and satisfying appropriate equation in 
R x (0, oo) and an initial condition. 

In order to deal with non-cooperative system, we shall assume that there are 
additional two monotone operators f ± , one lies above and another below / 
with the corresponding equations 

u t = Du xx + f + (u),xeB.,t>0. (2.5) 

tt( = D« n + f(u),i6l,i>0. (2.6) 

Such an assumption will enable us to make use of the corresponding results 
for cooperative systems in [T8|l3~l~] to establish spreading speeds for (12.31) . 

(HI) i. Let /, f ± : W N — > W N be Lipschitz continuous, twice piecewise contin- 
uous differentiable function such that 



r(u)<f(u)<f+(u),ue 



DiV 



ii. Let << k- = (fcf) < k = (h) < k + and /(0) = f(k) = and 
assume that there is no other positive equilibrium of / between and 



k (that is, there is no constant v ^ k such that f(v) = 0, << v < k). 
/ ± (0) = f ± {k ± ) = 0. There is no other positive equilibrium of / 
between and k . 
iii. ( 12. 5 p and ( 12. 6 p are cooperative (i.e. diff{u) > for u G [0, k + ],i ^ j). 
/ ± (u), f(u) have the same Jacobian matrix /'(0) at w = 0. 

A traveling wave solution u of (12. 3p is a solution of the form u = u(x + ct), u G 
C(R,R N ) . Substituting u(x,t) = u(x + ct) into (Q and letting f = x + ct, 
we obtain the wave equation 

W(£) - cu'(£) + /(w(0) = 0, e € R. (2.7) 

Now if we look for a solution of the form (iii) = (e x ^r)\) , A > 0, 77^ = (rj\) >> 
for the linearization of (12. 7p at the origin, we arrive at the following system 
equation 

diag(^A 2 - cA)r/ A + f\0) Vx = 

which can be rewritten as the following eigenvalue problem 

-A X 7] X = crix, (2.8) 

where 

A x = (a^) = diag(^A 2 ) + /'(0) 

The matrix /'(0) has nonnegative off diagonal elements. In fact, there is a 
constant a such that f'(0)+al has nonnegative entries, where I is the identity 
matrix. 

By reordering the coordinates, we can assume that /'(0) is in block lower 
triangular form, in which all the diagonal blocks are irreducible or 1 by 1 
zero matrix. A matrix is irreducible if it is not similar to a lower triangular 
block matrix with two blocks via a permutation. Let p(A) be the spectral 
radius of A. For a matrix A with nonnegative off diagonal elements, from the 
Perron- Frobenius theorem, we shall call the eigenvalue 

■qj(A) = p(A + aI) -a 

of A, which has the same eigenvector, the principal eigenvalue of A (see e.g. 
[TTf3Tj ). Here A + al is nonnegative, and p(A + al) is the spectral radius of 
A + al. 

We make assumptions on A x and also requires / grows less than its lineariza- 
tion along the particular function v x e~ Xx [31]. Such a condition can be satisfied 
for many biological systems. 

(H2) i Assume that A x is in block lower triangular form and the first diag- 
onal block has the positive principal eigenvalue ty(A\), and *£?(A X ) is 



strictly larger than the principal eigenvalues of all other diagonal blocks 
for some interval [0, A*] of A. Assume that for A G [0, A*] (we also as- 
sume that A c » < A* < oo, see Lemma [2.11 for A c *), there is a positive 
eigenvector v\ = (u\) » of A\ corresponding to ty(A\). 
ii Assume that for each A G [0,A*], 9 > 

^(K) < f'^Wx- 



Let 



$(A) = i*(A A ) > 0. 



According to Lemma I2.1[ we can expect the graph of $ as in Fig. LTJ For the 
example in Section IH $ is a strictly convex function of A and, clearly satisfies 
Lemma 12.11 



Now we state Lemma 12.11 which is a analogous result in Weinberger [3i 
and Lui [18J. However, due to the fact that /'(0) is only quasi-positive and 
the elements of A\ are not necessarily log convex, some of its proof here are 
different from Lui [18]. A similar result is included in [30]. A theorem on the 
convexity of the dominant eigenvalue of matrices due to Cohen [I] is used to 
show that ty(A\) is convex function of A. Lemma [2.11 improves [3"T| Theorem 
4.2] by eliminating the case (b) in |31[ Theorem 4.2]. 

Lemma 2.1 Assume that (HI) — (H2) hold. Then 

(1) $(A) -► oo as A ->■ 0; 

(2) If A* = oo ; $(A) ->> oo as A -»■ oo; 
^ $(A) zs decreasing as A = + ; 

(^ \I/(A0 zs a convex function of X > 0; 

f5j $'(A) changes sign at most once on (0,oo) 

^ $(A) /ias i/ie minimum 

c* = inf $(A) > 

A>0 V ' 

at a finite A c * > 0. 
(^7j For each c > c* , there exist a positive A c < A c * and 7 G (l,min{2, -f 1 }) 
such that 

$(A C ) = c, $(tA c ) < c. 



T/iai is 



and 



1 
t-A\ c ^a c = $(A c )z/ Ac = c^ Ac 



-v4 7Ac ^ 7Ac = $(7A c )z/ 7Ac < CI/. 



7 A C 

where v^ c » 0, v 7 a c >> are positive eigenvectors of j-A^ c ,-^-A^\ c 
corresponding to eigenvalues $(A C ) and $(7A C ) respectively. 



PROOF. The proof of the convexity of ty(A\) is similar to that in Crooks [5] 
for matrices with positive off-diagonal elements. It is easily seen that ty(A\) = 
p{A\ + al) — a is non- decreasing function of A > ([TTJ Theorem 8.1.18]). 
Further, a theorem on the convexity of the dominant eigenvalue of matrices 
due to Cohen jl] states that for any positive diagonal matrices Di,D 2 and 
* 6(0,1), 

*(*£>! + (1 + t)D 2 + /'(0)) < m(D 1 + f (0)) + (1 - t)^(D 2 + /'(0)) 

as before, here *&{A) is the principle eigenvalue of A. Now if a±,a 2 £ K and 
£6(0,1), 

(to x + (1 - t)a 2 ) 2 < ta\ + (1 - t)c^. 

This implies that 

*(A Ul+(1 _ t)A2 ) = *((tAi + (1 - £)A 2 ) 2 D + /'(0)) 
<M>(£A 2 D+(l-i)A 2 D + f(0)) 
<tM/(A 2 D + /(0)) 

+ (l-t)^(A 2 D + /'(0)) 
= t*(A Al ) + (l-t)*(A Aa ) 

Since ^(A A ) is a simple root of the characteristic equation of an irreducible 
block, it can be shown that *£?(A\) is twice continuously differentiable on K. 
Thus 

tf"(A) > 

and a calculation shows 

[A$(A)]' = #'(A) 

$'(A) = ~[*'(A)-$(A)] 

and 

(A 2 $'(A))' = A#"(A) > 0. 

As for (2), we need to prove that liniA_>oo A = oo. In fact, there exists an 
e > such that all diagonal elements of D — el are strictly positive, then 
^>(D — el) > and choose A large enough so that 

V(A X ) = V(D\ 2 + /'(0)) 

= v((D-eI) A 2 + (e\ 2 I + f 
> ^((L>-e/)A 2 
= A 2 *(.D-eI 

Thus limA^oo A = oo. (6) is a consequence of (l)-(5). (7) is a direct con- 
sequence of (l)-(6). It is just the fact that v\ » is a eigenvector of jA x 
corresponding to eigenvalue 3>(A\) for A = A c and 7A C . □ 




Fig. 1. The red curve is <1?(A). The minimum of $(A) is c* . For c > c* , the left 
solution of ^(A) = c is A c . 

We now recall results on the spreading speeds in Weinberger, Lewis and Li 
[3T] and Lui [IS]. While Theorem 4.1 [31] holds for non cooperative reaction- 
diffusion systems, it does require that the reaction-diffusion system has a sin- 
gle speed. In general, such a condition is very difficult to verify. In the same 
section, for cooperative systems, Theorem 4.2 in [31J provides sufficient con- 
ditions to have a single speed. The following theorem combines the results of 
Theorems 4.1 and 4.2 in [31], which can be a consequence of Theorems 3.1 
and 3.2 for discrete-time recursions in Lui 1181. 



Theorem 2.2 (Weinberger, Lewis and Li [31]) Assume (H1) — (H2) hold and 
A 2.3}) is cooperative. Then the following statements are valid: 

(i) For any uq G Ck with compact support, let u(x,t) be the solution of A 2. 3}) 
with {2.J$ . Then 

lim sup u(x, i) = 0, for c > c* 



t— >OD 



|x|>ct 



(ii) For any strictly positive vector u G WL N , there is a positive R^ with the 
property that if uq G Ck and uq > u on an interval of length 2_R W , then 
the solution u(x,t) of A2.3}) with \2.4\) satisfies 

liminf inf u(x, t) = k, for < c < c* 

tr+OC \ X \<tC 

In another paper [IB], for cooperative systems, Li, Weinberger and Lewis es- 
tablished that the slowest spreading speed c* can always be characterized 
as the slowest speed of a family of traveling waves. These results describe the 
properties of spreading speed c* for monotone systems. Based on these spread- 
ing results for cooperative systems, we will discuss analogous spreading speed 



results for non cooperative systems. 

3 Results on general non-cooperative systems 



In this section, we state a theorem for general partially cooperative reaction- 
diffusion systems ( 12 .3p which establishes the existence of traveling waves and 
spreading speed for a large class of non-cooperative systems. As we discussed in 
Section [21 assumptions (H1-H2) and the proofs in Section are different from 
those in [30J and the assumptions seem easier to verify. Although the existence 
of traveling wave solutions for cooperative systems are known (see, e.g. [16]), 
we shall prove the existence of traveling wave solutions for both cooperative 
and non cooperative systems as our proofs for non cooperative systems are 
based on those for cooperative systems. Further, in additions to the existence 
of traveling wave solutions, we shall be able to obtain asymptotic behavior of 
the traveling wave solutions in terms of eigenvalues and eigenvectors for both 
cooperative and non cooperative systems. The following theorem is our results 
for general non-cooperative reaction-diffusion systems. 

Theorem 3.1 Assume (HI) — (H2) hold. Then the following statements are 
valid: 

(i.) For any u G Ck with compact support, the solution u(x,t) of 112. 3\) with 
l\2.4\) satisfies 



lim sup u(x, t) = 0, for c > c* 

t ^ 00 \x\>tc 

(ii.) For any vector u G K , u » 0, there is a positive R u with the property 
that if uq G Ck and uq > u on an interval of length 2R U1 , then the solution 
u(x, t) of A2.3\) with \2.1$ satisfies 

k~ < liminf inf u(x, t) < k + , for < c < c* 

i->oo \x\<tc 

(Hi.) For each c > c* Ii2.3\) admits a traveling wave solution u = u(x + ct) such 
thatO ««(£) < fc + ,(el, 

AT < liminf w(£) < limsupu(£) < k + 

and 

lim u{i)e~ A ^ = v Kc . (3.9) 

If, in addition, Ii2.3\) is cooperative in Ck, then u is nondecreasing on R. 
(iv.) For c = c* /i2.y\) admits a nonconstant traveling wave solution u = u(x + 
ct) such that < u(£) <fc,(el, 

k~ < liminf u(£) < limsupn(^) < k + . 



(v.) ForO < c < c* Ii2.3\) does not admit a traveling wave solution u = u(x+ct) 
with liminf^oo -u(£) >> and w(— oo) = 0. 

Remark 3.2 In many cases, / can be taken as piecewise functions consisting 
of / and appropriate constants as demonstrated in Section HI In order to have 
a better estimate for the traveling wave solution u for non cooperative systems, 
it is desirable to choose two function f^ 1 which are close enough. The smallest 
monotone function above / and the largest monotone function below / are 
natural choices of f ± if they satisfy other requirements, See [2TIfl2|ll7] for the 
discussion for scalar cases and [55] for a partially cooperative reaction-diffusion 
system. Our construct of /~ in Section H] is different from the previous papers. 

Remark 3.3 When (12. 3p is cooperative in C/~, then f ± = f. 

Remark 3.4 Assumptions (Hl)(i-ii) imply that C^+ is an invariant set of 
(12. 3 p in the sense that for any given it G Ck+, the solution of (12 .3j) with the 
initial condition Uq exists and remains in Ck+ for t G [0, oo). In fact, for a 
given uq G C&+, let u(x, t) be the solution of (12. 3p with the initial condition 
uq. Theorem 15.11 implies that 

0< u(x,t) <fc + ,i6l,t>0. 

Now according to Smoller [251 Theorem 14.4] (12. 3 j) (and also (12. 5p . (J2.6J1 ) has 
a solution u for t G [0, oo) and < u < k + if the initial value uq is uniformly 
continuous on R. 

We shall prove Theorem 13.11 (i)-(iii) in Section |5] and (iii)-(v) in Section |6j 



4 Results on a model arising from epidermal wound healing 



In Section HI we shall apply the general results in Section [3] to the model (II. ip 
arising from epidermal wound healing. This model is not cooperative because 
of the fact that h(u\) is not monotone. We shall establish the existence of 
traveling waves as well as the results on the speed of propagation to (II. ip . 
In addition, we characterize the spreading speed as the slowest speed of a 
family of non-constant traveling wave solutions of (II. ip . The spreading speed 
for (II .ip was discussed in [24~|IT9] based on numerical methods and singular 
perturbation techniques for several special cases, for example, d\ = 0. 

Recall that h(ui) = ~^tA ■ It is easy to (II. ip has two equilibria (0,0) and 



10 



[1,1). In fact, the following equalities hold at its non-trivial equilibrium 

2^ + K- 1 

Kui) = — ^ (410) 

u 2 = h{ui). 



Now it is clear that ( I4.10p has only one positive solution (1,1). In fact, 

in (14.10)) is increasing and convex on (0, oo) and the first equation of (14.1 Op 

has only one solution U\ — l. 

We now need to check (H2). The linearization of ( II. ip at the origin is 



diAui + (1 — 2k)Ui 
at 

% = d 2 Au 2 + b(h'(0) Ul - u 2 ) 
at 

where h'(0) = ±#. The matrix A x in (ESD for flU)) is 



(4.11) 



C 2 

/ diA 2 + l-2K 



A, = (a^) 



\^'(0) d 2 X 2 -b 



(4.12) 



It is easy to see that 

h(ui)<ti(0)ui, uie(0,oo). (4.13) 

In order to use Theorem 13.11 we shall define the two monotone systems. Note 
that h{u\) achieves its maximum value -^$- when u\ = £. 

h(ui), < u\ < C, 

' h((), Ul > c. 

and the corresponding cooperative system is 



diAui + s{u 2 )ux{2 - Mi) - ui 
at 

d 2 Au 2 + 6(/i + (ui) - u 2 ) 



(4.14) 



dt 
In a similar manner, one can find (I4.14p has two equilibrium (0, 0) and (kf, k 2 ) 

11 




satisfying 



Fig. 2. The construction of h + and h . The red curve is h. 



h + (kt) = ^ 



(4.15) 



k+ = h + (kt 



Since h + = h for u\ < £, then A;^ > £ (if kf < ( , then kf = 1 and 
C > 1) and h + (ki) = -^-. Solving fcj*" directly from f |4.15|) gives that A;^ = 

> 1 > C- It follows that kt = h + (kt) = i±£ > 1. 



i+(^-i)" 



2C 



Now there is a /i G (0,£] such that h(h ) < mm{l, h(kf)} and define 



h (ui) 



h(ui), < u\ < ho, 
h(ho), mi > ho- 



lt is clear that 

< h~{u x ) < h(u x ) < h + ( Ul ) < h'^juuux G (0,fc+] 
and h~{u\) < 1 for U\ > 0. 
The corresponding cooperative system for h~ is 
du i 



du 2 



diAui + s(u 2 )ui(2 — ui) — U\ 
d 2 Au 2 + b(h~(ui) - u 2 ) 



(4.16) 



In a similar manner, one can find (14. 16f) has two equilibrium (0, 0) and (k ± , k 2 



12 



satisfying 

h ( fc i ) = — — K (4.17) 

k 2 = h~(fo). 

Since h~ = h for u\ < ho, then k^ > h and h~(ki) = h(h ). Solving k± 
directly from (J4.17P gives that fcf = 2 — 1+(fe / fe 1 n_ 1 n k < 1 as /i(/io) < 1- On the 
other hand, because h(ho) > 0, a simple calculation shows that k{ > y^7 > 0. 
As before we have < k 2 = h~(ki) = h(ho) < 1. 

Thus, 

(0,0) << (/cf;,^) < (1,1) < (/c+,A:+). 

We can always extend (l4.14IUI4.16p to be Liptschz continuous in M 2 without 
changing the functions in the region [0, k + ]. Then Theorem 15.11 implies that 

<u(x,t) <k + ,xeR,t>0. 

Thus we are only interested in the invariant region. Now it is straightforward 
to check all other conditions of (HI) (i)-(iii) . 

The spreading results for the cooperative systems were used to establish in [35] . 
We now demonstrate Theorem 13.11 can be used to establish spreading speed 
and traveling wave solutions of the nonmonotone system fll.lf) and summarize 
our main results in the following Theorem. 

Theorem 4.1 Let d\, d 2 be all positive numbers and k G (0, |), £ 6 (0, 1), 

do b , 

i <2+ —. < 418 > 

and 

( . 1-2/c . 

2«tf (0) l + —f— dl ^ rf2 ' , , 

, V ' <{ a i _V (4.19) 

"}(2-|)^ + l,rfi<rf 2 
Then the conclusions of Theorem \3. 1\ hold for U.l\) where the minimum speed 

c- v /c 2 -4di(l-2/c) 



c* = 2^/(1 - 2k) di, A c = — * — 2dT^ > ° anc ^ ^ ^ s defined in \4-21\) . 



That is, the solution (ui(x,t),U2{x,t) of ( ti. 1\) satisfies 

(i.) If the functions (ui(x, 0), U2(x, 0)) < (fci,/c 2 ) are nonnegative continuous 
and have compact support, then 

lim sup (ui(x, £), «2 (#,£)) = (0, 0) for c> c* 
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(ii.) If the functions (u\(x, 0),U2(x,0)) < (^1,^2) are nonnegative continuous 
and Ui(x, 0) ^ ; then 

(k± , k 2 ) < liminf inf (ui(x, t),Wx, t)) < (ki, k^), for0<c<c* 

t->-oo \ x \<tc 

(Hi.) For each c > c* M.l\) admits a traveling wave solution («i(£) , w 2 (£)) sitc/i 
&a* (0,0) « («i(0,«2(0) < (#,#),££ R, 

(K,K) < liminf(«i(£),«2(£)) < limsup(tti(£),U2(0) ^ (kf^t) 

and 

Km K(0,u 2 (0)e- Ac? = ^A c . (4.20) 

4— >— 00 

(w.j Fore = c* ( fi?.ff|) admits a nonconstant traveling wave solution (ui(£), n 2 (£)) 
such that (0,0) « («i(f),«2(0) ^ (^,^).^1, 

(^r,^) < liminf(«i(0,tt2(0) ^ lim sup «(^) < (k+,k£). 

(v.j For < c < c* ( QO[) does not admit a traveling wave solution (wi(£), w 2 (£)) 
to£/j liminf£_ KX) ('Ui(£), w 2 (£)) >> (0)0) anc ^ ( u i( — °°)> u 2(— 00)) = 0. 

Now we shall verify (H2) for (j4.14p . In fact, the principle eigenvalue of A\ is 
ty(A\) = di\ 2 + 1 — 2k, which is a convex function of A. Furthermore, 

$( A ) = E^) = ^ A2 + l - 2k 

A A 

satisfies the results of Lemma I2TT1 In fact $(A) is also a strictly convex function 



of A. The minimum of $(A) is c* = 2i/(l - 2k) di when A = ^^ K . For each 
c > c*, the left positive solution of $(A) = c in Lemma [2.11 is 



c-^c 2 -4di(l-2«) 

Ac= 2d^ 

In particular, 



A c * = - — 1^—. 
Va\ 

For each < A < A c », the positive eigenvector of A\ corresponding to ^(A\) 
is 

Because of ( 14.181) . ^ is a strictly positive vector for A G [0, A c *]. This is clear 
when di > d 2 . If di < d2 and (14.181) holds, for < A < A c . , we have 

(di - d 2 )A 2 + 1 - 2k + 6 > (di - d 2 )A^ + 1 - 2k + b > 0. 
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Further from (j4.2ip we can see that 

v\_ bh'(0) fe'(Q) 

v\ ~ (dx - d 2 )X 2 + 1 - 2k + b ~ a 

(d 1 -d 2 )\ 2 +l-2K t? x r n a 1 ., - , ,, , 1 , 1-2k 



where a = 1 + igizg^+iz^ . For A G [0,A C .], it is clear that a > 1 + ±=p if 

d 2 N 1-2, 
di' ft 



di > d 2 and ct > (2 - #)^ + 1 > if d x < d 2 and (Qgj) holds. 



Let 



( Wl ,u 2 ) = (0,0— )»(0,0), #>0. 



Thus (H2)(ii) for (j4.14p is equivalent to the following two inequalities 

(ku 2 + 1 - k)ui(2 - Mi) - tti < (1 - 2«)ui , . 

bh(u 1 )-bu 2 <bh'(0)u 1 -bu 2 . 

Because h(ui) < h'(0)ui, (j4.22p is equivalent to the following inequality 

(2-ui)k<(1-k)— (4.23) 

u 2 

and the following inequality suffices to verify ( I4.23P 

2k a 



h'(0) 



which is true with ( I4.19P because of the estimates for a for d± > d 2 and 
d\ < d 2 . Notice that h^ and h are identical around the origin. By the exact 
same arguments (just replacing kf by k~ and h + by h~), we can verify that 
(H2) holds for (OEjl as well. 

It remains to show that the condition (ii) in Theorem 13.11 can be satisfied if 
u±(x, 0) ^ 0. The arguments here is the same as in Weinberger, Kawasaki and 
Shigesada [35] . We choose positive constants p, r] so small that 

- d x p 2 + (1 - 2k) - (1 - k)t] > 0. (4.24) 

Since u 2 (x, t) > 0, we have 

— — > diAui + (1 - k)ui(2 -u\)—u\ 
at (4.25) 

= diAu\ + «i((l — 2k) — (1 — k)uij 

By the strong maximum principle we have u\(x, t) > for t > 0. Thus we can 
require that r) < Ui(x,ti) for some t\ > and |x| < j- and some t\ > by 
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Fig. 3. The simulations of the traveling wave solutions of (jl.ip . We choose 
d 1 = d 2 = l,K = 0.05, C = 0.1, b = 0.05. 



choosing r/ small enough. If {u\, u 2 ) is the solution of 



dui 
du 2 



diAu\ + ui((l — 2k) — (1 — k)u\ 



d 2 Au 2 + b(h (ui) - u 2 ) 



(4.26) 



with the initial values 



uAxM) 



rjcos(px) for |x| < ^ 




2/0 



for \x\ > 7T 



n 2 (a;,ti 



(4.27) 



It is clear that (I4.26P has two equilibriums (0,0) and (^f,^(^f)) >> 0. 
Furthermore, there is no other stationary solution of (I4.26P between the two 
equilibriums. 



The comparison principle shows that the components of (ui,u 2 ) are lower 
bounds for (m, u 2 ) when t>t x . The inequality (14.241) shows that both (^-, ^) 
are nonnegative at t — t\, and the comparison principle then implies that 
(ui,u 2 ) are nondecreasing in t. It follows that (ui,u 2 ) monotonically con- 
verges to (-j^r 1 , /i~( 1 1 ~ 2w )) uniformly in x on every bounded x- interval. Because 
(ui,u 2 ) > (ui,u 2 ), it follows that if we choose two positive constants (cji,Ui) 



_ . ,. , 1 _^ 1 )) then for all sufficiently large t, the condition 
(ii) in Theorem 13.11 is automatically satisfied on the fixed interval |x| < 1R^. 
We thus obtain the statement (4.21) without an extra condition. 



Fig. [3] are the simulations of the traveling wave solutions of (II. ip . We choose 
d\ = d 2 = 1,k = 0.05, ( = 0.1, b = 0.05, which satisfy the conditions of 
Theorem 14.11 Note that the traveling solutions are not monotone and the 
minimum speed c* = 2-i/(l — 2K)d± = 1.89. 
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5 The Spreading Speed 

5.1 Comparison Principle 



We state the following comparison theorem for reaction-diffusion equations in 
File |H] . The comparison principle is a consequence of the maximum principle 
(see, e.g., Protter and Weingberger [13]). Another form of the comparison 
principle from [35] was used in [30]. As a result, the proof in Section [5] is 
somewhat different from [30]. By a suitable modification of the functions, the 
conditions in this paper seem easier to verify. 

Theorem 5.1 Let D be a positive definite diagonal matrix and F = (Fj) 
vector-valued continuous functions in M. N . Assume that i 7 ^ = (F- ) are Lips- 
chitz continuous on R and F^ 1 is nondecreasing in all but the jth component, 
j = 1, ..., N, and satisfying 



F + {u) > F{u) > F~(u),ue 



)iV 



Let u^u^ be continuous functions from R x (0,T) into M, N , C 2 in R x (0,T), 
bounded and satisfying, for i — 1, ..., n 

u t - Du xx = F(u) 

u+-Du + xx >F + {u+) (5.28) 

Ut -Du~ x < F~(u~) 

Ifu + (x,to) > u(x,to) > u~(x,to), i6l, then 

u + (x,t) > u(x,t) > u~(x,t), x E R,£ G [t ,T]. 

We are now able to prove Parts (i) and (ii) of Theorem 13.11 

5.2 Proof of Parts (i) and (ii) of Theorem \3.1\ 

Part (i). For a given uq G Ck with compact support, let u + (x, t) be the solutions 
of (12.51) with the same initial condition w , then Theorem 15.11 implies that 
u + (x, t) G Ck+ and 

< u(x, t) < u + (x, t), x G R, t > 0. 

Thus for any c > c*, it follows from Theorem 12.21 (i) that 

lim sup u + (x,t) = 0, 

*"*' 00 \x\>tc 
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and hence 



lim sup u(x, t) = 0, 



*^°° \x\>tc 



Part (ii). According to Theorem 12.21 (ii). for any strictly positive constant u>, 
there is a positive R u (choose the larger one between the R u for (I2.5P and the 
Ruj for (12. 6p ) with the property that if u > u on an interval of length 2R U]1 
then the solutions w ± (a;, t) of (I2.5P and (12. 6p with the same initial value u® are 
in Ck+ and satisfy 

liminf inf u (x,t) = k , for < c < c*. 

t^oo \ x \<tc 

Thus, Theorem 15.11 implies that 

u~(x,t) < u(x,t) < u + (x,t),x eR,t>0. 
Thus for any c < c*, it follow from Theorem 12.21 (ii) that 

liminf inf u (x,t) = k , 

t^oo \ x \<ct 

and hence 

k~ < liminf inf u(x,t) < k + . 

t-S>oo \ x \<ct 

a 



6 The characterization of c* as the slowest speeds of traveling waves 

6. 1 Equivalent integral equations and their upper and lower solutions 

In order to establish the existence of travel wave solutions, we fist set up 
equivalent integral equations. Similar equivalent integral equations were also 
used before, see for example, Wu and Zou [37], Ma [2~Tf22] and the author [28] . 
For the convenience of analysis, in this paper and [2S], both Aij, A 2 j are chosen 
to be positive, and — A^, A 2 j are solutions of (16. 29 p . 

Let (3 > max{\difj(x)\,x G [0,k + ],i,j = l,...,iV} > 0. For c > c*, the two 
solutions of the following equations, 

di\ 2 -cX- /3 = 0,i = l,...,N (6.29) 

are —\u and A 2 « where 

Xu ~ 2d > °' A2J " U > °- 



(6.31) 



We choose /3 sufficiently large so that 

X 2i >X u >2A c ,i = l,...,N. (6.30) 

Let u = (ui) G Ck and define a operator T[u] = (7i[u\) by 

a i\*li "T ^2i) Woo 

+ f°° e M ^- s) H t {u{s))ds\ 

where 

H l (u(s))=(3u l (s) + f l (u(s)), 
Ti[u],i = 1,...,N is defined on R if Hi(u),i = 1,2 is a bounded continuous 
function. In fact, the following identity holds 

e- Xl ' { ^ s) (3ds+ / e X2 ^-*)pds 



di(Xu + X21) W-00 Jt, 

P / 1 , 1 \ = I (6.32) 

<MAii + A2i) ^Aij A2r <^i(AijA2j) 
= 1. 

In fact, a fixed point w of T or solution of the equation 

u(0=T[tt](0 ^Gl, (6.33) 

is a traveling wave solution of (12. 3p in Lemma 16.11 Lemma 16.11 summarizes 
the conclusion, which can be verified in the same manner as in [28] for scalar 
cases by directly substituting derivatives of u(£) into (12. 7p . We omit its proof 
here. 

Lemma 6.1 Assume (HI — H2) hold. If u G Ck is a fixed point ofT[u], 

w(£)=T[w](0 teR, 
then u is a solution of \2.Tty . 

We now define upper and lower solutions of (I6.33J) . <p + and <p~ , which are only 
continuous on R. Similar upper and lower solutions have been frequently used 
in the literatures. See Diekmann [7], Weinberger [31], Liu [18], Weinberger, 
Lewis and Li [21], Rass and Radcliffe [20], Weng and Zhao [SB] and more 
recently, Ma [22], Fang and Zhao [8] and Wang J2SH2H]- In particular, it is 
believed that the vector- valued lower solutions of the form in this paper first 
appeared in [36J for multi-type SIS epidemic models. In this paper, the upper 
and lower solutions here are defined and verified for the integral operator other 
than differential equations. 
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Fig. 4. For each i, the curve above is (fif and the below is (\> i . 

Definition 6.2 A bounded continuous function u = (m) G C(R, [0, oo)^) is 
an upper solution of (16.331) if 

Ti[u](0 < m(Z), for aU £ G R,i = 1, ..., AT; 

a bounded continuous function u = (uj) G C(R, [0, oo)^) is a lower solution 
of (RT35JI if 

Ti[u}(0 > Ui (0, for all £ G R,i = 1,...,JV. 

Let c > c* and consider the positive eigenvalue A c and corresponding eigen- 
vector v^ c = (u l A ), 7 in Lemma [2.11 and q > 1. Define 

+ (O = (</>+), 

where 

0+ = min^, ^ c e Ac «}, i = 1, ..., N, £ G R; 

and 

4>-(0 = (<PT), 

fa = max{0, z/l c e Ac « - ^e^}, i = 1, ..., A, £ G R. 

It is clear that if £ > -^ 0+(£) = A*, and £ < -j^ 6 -, <#(0 = A^,i = 
1,...,N. 

Similarly, if C > (1 _ 7 fc e > 0< (0 = °> and for £ < ( i- 7 fc e , 

0. r (O = ^ Ac€ - g^e^.i = 1, ...,iV. 
We choose g > 1 large enough that 

ln(g^) IhjJl 
*£ < *£. i — l AT 

(1-7)A C A c ' '"■' 

and then 



<Pt(0><P7(0^ = h-,N,teR. 
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We now prove that 4> + and 4> are upper and lower solution of (I6.33|) respec- 
tively. 

Lemma 6.3 Assume (HI) — (H2) hold and 112. fy) is cooperative. For any c > 
c* , (p + defined above is an upper solution of Ii6.33\) . 



PROOF. Note that when (J2.3P is cooperative, f = f (see Remark [3731) . Let 

£* = A Ac , i = 1, ..., N. We first need to verify that u = (ui) = <p + is an upper 
solution of (12.71) in the sense that 



di<(o - «4(o + fMO) < o, ,i = i, ..., i\r, ^ er- (6.34) 

In view of the fact that /*(«(£)) < /i(fc) = for £ > 4*, it is clear that (16.341) 
holds for £ > £*. For £ < 4*, ^i(0 = ^i c eAc? - From (H1-H2), we have, for 
f < 4* andi = 1,...,JV 

n 

i=i 
= E«K^-cA c ^ c e^ (6 ' 35) 

= A c $(A c )z/l c e Ac « - cA c ^ c e Ac « 
= 0. 

Since u'^*—) > u' i (C,*+) = 0, the proof of Lemma 2.5 in Ma [21] immediately 
implies that 

7TM£) < «*(£), for all £ 6 R,i = 1, ..., JV. 

Thus -u = (uj) = 4> + is an upper solution of (16.331) . Note that the upper 
solutions in this paper are defined for equivalent integral equations other than 
differential equations in [2TJ. □ 

We now need the following estimate on /, which is an application of the 
Taylor's Theorem for multi- variable functions. Also see |28j . 

Lemma 6.4 Assume (H1 — H2) hold. There exist positive constants bij,i,j = 
1, ...,N such that 

N N 

fi(u) > ^2djfi(0)uj - ^bij^Uj) 2 , u = (ui),Ui e [0,k],i = 1,...,N. 

i=i 3=1 

Lemma 6.5 Assume (HI) — (H2) hold. For any c > c* , <f>~ defined above is 
a lower solution of A6.33\) if q (which is independent of £) is sufficiently large. 
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PROOF. Let Let 4* = (1 _^ e ,i = l,...,N. We first need to verify that 
u = (ui) = (f)~ is a lower solution of (12.71) in the sense that 

dtu'ltt) - cafe) + fMO) > 0, ^ £,z = 1, ...,N. (6.36) 

It is clear that (jOgJ holds for f > g. For £ < £*, 

From Lemma [6.41 we have, for £ < £* and z = 1, ...,N 
d iU '!{£) - au^) + f t (u(0) 

n 

> d t Ay Ac e A ^ - ck c v\ c e^ + £ a 3 / 4 (0)<e A ^ 



J c 7Ac? 



- g(^( 7 A c )X Ac e^ - c 7 A c ^ Ac e^« + £ ^(0)^ Ac e 

- Me 2Ac? 



(6.37) 



where M, = Ej=i M^ij 2 > °- Note that e (2 - 7)Ac « is bounded above for 
£ < &• With the same argument in (16.351) . in particular, by Lemma 12.1 [ we 
get for sufficient large q and £ < & 

> g(c 7 A c - 7 A c $(7A c ))^ Ac e 7Ac? - M*e 2Ac « 

= g(c - $( 7 A c )) 7 A c z/; Ac e^ - M t e 2A ^ (6 . 38) 

> (q(c - $( 7 A c )) 7 A c ^ Ac - Mie^^y^* 
>0. 

Since i4(£* — ) < «$(£*+) = 0, the proof of Lemma 2.6 in Ma [21] immediately 
implies that 

Ti[u](0>Ui(0i foraU£eR,i = l,...,JV. 
Thus m = (wj) = + is a lower solution of (I6.33p . Note that the upper so- 
lutions in this paper are defined for equivalent integral equations other than 
differential equations as in [21]. □ 



6.2 Proof of Theorem \3.1\ (Hi) when Ii2.3\) is cooperative 



In this section, we assume that (12.31) is cooperative and prove Theorem 13.11 
(iii). As we note in Remark 13.31 / = /. Many results in this section are 
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standard. See, for example, Ma [21f22] and Wang [28|30] . Define the following 
Banach space 

S e = {u = ( Ui ) : ui G C(R),smp\ui(^)\e~^ < oo,z = 1,...,JV} 
equipped with weighted norm 

N 



U 



Q 



Esup|«i(0|e- rf , 



i=\ 



where C(R) is the set of all continuous functions on K and g is a positive 
constant such that g < A c . It follows that </) + G £ e and <p~ E £ g . Consider the 
following set 

A={u = ( Ul ) : Ui e C(R) e £ e AT(0 < ^ < #(£U e 1M = 1,...,JV.} 

We shall state the following lemmas. It is standard procedures to prove that 
T[u] is monotone, continuous and maps a bounded set in A into a compact set. 
We omit their proofs here. The proof of compactness can be carried out two 
steps. First we show that T is equicontinuous, and then Ascoli's theorem and 
standard diagonal process can be used to prove that T is relatively compact 
E e . The proofs of Lemmas I6.6|, 16.7^ I6.8l a.re almost identical to those in 



Lemma 6.6 Assume (HI) — (H2) hold and difj >Q,i=£j on [0, k\. Then T 
defined in A6.31\) is monotone and therefore T(A) Q A. Furthermore, Ti[u} is 
nondecreasing if u G A and all of Ui are nondecreasing. 

Lemma 6.7 Assume (HI) — (H2) hold. Then T ' : A^r £ e is continuous with 
the weighted norm ||.|| e . 

Lemma 6.8 Assume (HI) — (H2) hold. Then the set T(A) is relatively com- 
pact in £ e . 

Now we are in a position to prove Theorem 13.11 when (12. 3 j) is cooperative. 
Define the following iteration 



u 



1 = (ul) = T[0 + ], u n+1 = «) = T[u%n> 1. (6.39) 



From Lemmas 16. 3[ 16. 5[ 16. 6[ u n is nondecreasing on M and 

</>-(£) < < +1 (0 < K(0 < ttiZU e R, n > l,n = 1, ..., N. 

By Lemma [6781 and monotonicity of (u n ), there is u G A such that lim n _ s . 00 || , u n — 
u\\ e = 0. Lemma [6.71 implies that T[u] = u. Furthermore, u is nondecreasing. 
It is clear that \m\^_ 00 u i (^) = 0,i = 1,...,N. Assume that lim^oo Ui(£) = 
k'; L ,i = 1, ..., N k[ > 0, i — 1, ..., N because of u G A. Applying the Dominated 
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convergence theorem to (16.31 j) . we get k[ = \{fik\ + fi(k[, ■■■,k' n ) By (HI), 
k[ = ki. Finally, note that 

"l(e Ac( ~ qe^ < «,(£) < ^^G R. 
We immediately obtain 

lim Ui (Oe- A ^ = vl,t = l,...,N. (6.40) 

£— >— oo 

This completes the proof of Theorem 13.11 (iii) when (12 .3ft is cooperative. 
6. 3 Proof of Theorem \3.1\ (iii) 



PROOF. Theorem 13.11 (iii) is proved when (12. 3 p is cooperative in the last 
section. Now we need to prove it in the general case (12 .3p is not necessar- 
ily cooperative. In order to find traveling waves for (12. 31) . we will apply the 
Schauder's fixed point theorem. Similar arguments can be found in the refer- 
ences, e.g. (22] and [30J. 

Let u = (ui) G A and define two integral operators 

7*[«] = (7f[u]) 
for /~ and / + 

7? [«](*) 



d>i(\u + A2J 
and 



e- Xl ^s)H±(u(s))ds+ / e X2 ^- s) H^{u{s))ds] 

-00 J£ 



(6.41) 



Ht(u(s)) = Pu t (s) + ft(u { s)). 

As in Section I6T2"} both T + and T~ are monotone. In view of Section I6T21 and 
the fact that /~ is nondecreasing, there exists a nondecreasing fixed point 
u~ = (u^) of T~ such that T~[u~] = vr , lim^oo u~{£) = k~,i = 1,...,N, 
and lim^^oo u^(£) = 0,i = 1,...,N. Furthermore, \im^_ 00 u~(^)e~ Ac ^ = 
v\ , i = 1, ..., N. According to Lemma 16. 3\ <p + (with k being replaced by k^) 
is also a upper solution of T ± because the proof of Lemma 16.31 is still valid if 
/ is replaced by f ± . Let 

^(0 = (tf(0), 

where 

0+(£) = mm{kt, vle A ^} } i = l,...,N,£e R; 



24 



It follows that Ui(g) < fit, £ G R, i = 1, ..., N. Now let 

B = {u : u = (m) e £ e ,uT(£) < mi® < ${£),£ e (-00,00),*= 1, ...,N}, 

(6.42) 
where £ e is defined in Section 16.21 It is clear that B is a bounded nonempty 
closed convex subset in £ e . Furthermore, we have, for any u = (ui) G B 



U: 



%-[u-} < TT[u] < Ti[u] < % + [u] < 7- + [0+] < <Pt,i = 1, ...,N. 



Therefore, T : B — )■ 0. Note that the proofs of Lemmas 16. 7\ 16.81 are valid 
if (12.31) is not cooperative. In the same way as in Lemmas 16.71 16.81 , we can 
show that T : B — > B is continuous and maps bounded sets into compact sets. 
Therefore, the Schauder Fixed Point Theorem shows that the operator T has 
a fixed point u in B, which is a traveling wave solution of (12. 3p for c > c*. 

Since uj~(£) < Wi(£) < <f>t{Q> £ e ( — °°) °°); * — 1) ■••> -^ it is easy to see that 
for i = 1, ..., iV, lim^.oo u t -(£) = 0, lim^^ Ui(£)e 



-A c £ 



'A c > 



fc < liminf u(£) < limsup-u(£) < /c H 
and < u~(£) < «<(£) < fe+,£ G (-00,00). □ 



(5.^ Proof of Theorem \3.1\ (iv) 



PROOF. We adopt the limiting approach in [5] to prove Theorem 13.11 (iv). 
For each n G N, choose c n > c* such that lim^oo c n = c*. According to 
Theorem 13.11 (iii), for each c n there is a traveling wave solution u n = (uf) of 
(12T3J) such that 

and 

kr < liminf<(£) < limsup<(£) < kf,i = 1,...,N. 

It can be shown that (u n ) is equicontinuous and uniformly bounded on R (see, 
e.g. [2H Lemma 5.3]), the Ascoli's theorem implies that there is vector valued 
continuous function u = (u l ) on R and subsequence (w„ m ) of (u n ) such that 

Jim o u flm (0=«(0 

uniformly in £ on any compact interval of R. Further in view of the dominated 
convergence theorem we have 

u = T[u](fl. 
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Here the underlying Ah, A2J of T is dependent on c and continuous functions 
of c. Thus u is a traveling solution of (J2.3P for c = c*. Since, for each c n , u n G i3 
where S is defined in (16.421) . it is easy to see that u satisfies 

k~ < liminf w(£) < limsupw(^) < k + 

Because of the translation invariance of u n , we always can assume that u n (0) < 
~fc~ for all n. Consequently u is not a constant traveling solution of (12. 3p . □ 



6. 5 Proof of Theorem \3.1\ (v) 



PROOF. Suppose, by contradiction, that for some c G (0, c*), ( 12 .3p has a 
traveling wave u(x,t) = u(x + ct) with liminf^oo «(£) >> and u{— 00) = 0. 
Thus m(x, t) = u(x + ct) can be larger than a positive vector with arbitrary 
length. It follows from Theorem 13.11 (ii) 

liminf inf u(x,t) > k~ » 0, for < c < c* 

t^oo \ x \<ct V ' ~ 

Let c G (c, c*) and x = ct. Then 



lim u[ — (c — c)t) = lim u(—ct,t) > liminf inf u(x,t) » 0. 

However, 

lim it( — (c — c)t) = u(—oo) = 0, 

which is a contradiction. □ 
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